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Abstract 

We derive the contributions to the quantum transport equations for electroweak baryogenesis 
due to decays and inverse decays induced by tri-scalar and Yukawa interactions. In the Mini- 
mal Supersymmetric Standard Model (MSSM), these contributions give rise to couplings between 
Higgs and fermion supermultiplet densities, thereby communicating the effects of CP- violation in 
the Higgs sector to the baryon sector. We show that the decay and inverse decay-induced contri- 
butions that arise at zeroth order in the strong coupling, a s , can be substantially larger than the 
0(a s ) terms that are generated by scattering processes and that are usually assumed to dominate. 
We revisit the often-used approximation of fast Yukawa-induced processes and show that for re- 
alistic parameter choices it is not justified. We solve the resulting quantum transport equations 
numerically with special attention to the impact of Yukawa rates and study the dependence of the 
baryon-to-entropy ratio Yb on MSSM parameters. 
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I. INTRODUCTION 



The origin of the baryon asymmetry of the universe (BAU) remains an open question 
for particle physics, nuclear physics, and cosmology. Although the size of the BAU cannot 
be explained within the framework of the Standard Model (SM), there exist a variety of 
SM extensions that may allow for successful baryogenesis. Scenarios in which the BAU 
is produced at the electroweak phase transition are particularly attractive since they can 
be tested with laboratory experiments. To the extent that the masses of the particles 
responsible for baryogenesis are not too different from the weak scale, their dynamics can be 
studied using a combination of collider experiments, precision electroweak measurements, 
and CP-violation studies. 

In order to carry out robust tests of electroweak baryogenesis (EWB), it is necessary 
to delineate systematically the quantitative relationship between EWB and experimentally 
accessible observables. The motivation for doing so has been heightened by the prospect 
of significant new experimental information in the near term. Studies at the Large Hadron 
Collider will search for the existence of new particles at the TeV scale. At the same time, a 
new generation of searches for permanent electric dipole moments of the electron, neutron, 
and neutral atoms will look for the effects of "new" CP-violation with several orders of 
magnitude better sensitivity than given by current experimental limits (see, e.g., Refs. 0,0] 
and references therein). Should either the LHC or EDM searches discover evidence for new 
physics at the electroweak scale, then precision studies at both the International Linear 
Collider and low-energy facilities should provide detailed information about the structure 
of the new physics. To the extent that the theoretical treatment of EWB is on sufficiently 
firm ground, these experimental efforts may either confirm or rule out this paradigm for the 
BAU. 

The basic physical picture of EWB was developed over a decade ago 0, 0, II 0] ( see 
for a review). The elements include a first-order electroweak phase transition, in which 
bubbles of broken electroweak symmetry expand and fill the universe as it cools through the 
transition temperature. CP- and C-violating interactions between fields in the plasma at 
the phase boundary create a net chiral charge that is injected into the region of unbroken 
electroweak symmetry, driving the weak sphaleron processes that create non-zero baryon 
number density, rig. The expanding bubbles then capture the non-zero % in the region 
of broken electroweak symmetry, where weak sphaleron processes are highly suppressed 
and unable to affect ub appreciably. It is crucial that the first-order phase transition be 
sufficiently strong in order to preclude "wash out" of non-zero baryon number. 

The earliest analyses based on this picture employed conventional transport theory to 
compute the production, diffusion, and relaxation of chiral charge at the phase boundary. 
Several groups have subsequently endeavored to put these computations on a more sophis- 
ticated footingbv using non-equilibrium quantum field theory techniques. As first pointed 
out by Riotto jalSl; only a non-equilibrium field-theoretic formulation can properly account 
for the quantum nature of CP violation as well as the decoherence effects due to the pres- 
ence of spacetime- varying background fields and the thermal bath of particles at the phase 
boundary. Using these methods Riotto 0] observed that conventional treatments may over- 
look significant enhancement of the CP-violating source terms in the transport equations 
associated with memory effects in the plasma. The presence of such enhancements could 
relax the requirements on new CP-violation needed for successful EWB, thereby allowing 
for consistency between the BAU and considerably smaller EDMs than previously thought. 
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This work was followed by the authors of Refs. 111! ], who adopted a similar approach 
to that of Ref. 0] in computing the CP-violating source terms while carrying out a more 
comprehensive phenomenology. The analyses of both groups were performed within the 
Minimal Supersymmetric Standard Model (MSSM). The non-equilibrium approach has also 
been pursued in Refs. 12[ HI 3 • 



Recently, we investigated the CP-conserving terms as well as the CP-violating sources 
in the transport equations using non-equilibrium field theory methods 0]. We found that 
there exists a hierarchy of physical scales associated with the electroweak phase transition 
dynamics that allows one to derive the transport equations from the Closed Time Path 
Schwinger-Dyson equations using a systematic expansion in scale ratios. Again in the MSSM, 
we computed the CP-violating sources and leading CP-conserving chiral relaxation terms 
associated with interactions of fermion and Higgs superfields with the spacetime-varying 
Higgs vacuum expectation values ( vevs) . Our results for the sources were consistent with 
those obtained in previous work 0,0,0, 3, 17], but we also found that enhancements in 



the relaxation rates could mitigate the effect of enhancements in the sources. 

A number of other contributions to the transport equations remain to be analyzed using 
non-equilibrium methods. Here, we focus on terms that link the dynamics of the quark 
supermultiplets with those of the Higgs scalars and their Higgsino superpartners. Impor- 
tantly, these terms are responsible for communicating CP-violating effects in the Higgs su- 
permultiplet densities to the quark supermultiplet densities, thereby allowing CP-violating 
interactions in the Higgs sector to contribute to baryogenesis. In the MSSM, the require- 
ment of a strong first-order phase transition (shown in [l8| to occur in the presence of a 
light right-handed stop) and constraints from precision electroweak data (requiring the left- 
handed stop to be heavy 0|) imply that it is the CP- violating interactions of the Higgs 
superfields - rather than those directly involving the squarks - that drives baryogenesis via 
this coupling between the two sectors. In extensions of the MSSM, such as the NMSSM 
or U(l)' models, the phenomenological requirements that preclude large effects from CP- 
violation in the squark sector can be relaxed and in this case it is important to know 
the relative importance of Higgs sector CP-violation. In either case, an analysis of the dy- 
namics whereby the baryon and Higgs sectors communicate is an important component of 
a systematic, quantitative treatment of EWB. 

Before providing the details of our study, we summarize the primary results, using the 
transport equation for the Higgs + Higgsino densities for illustration: 

r H . = - r B r JQT + B\-t r (° * + . (1) 



^ k H \k Q k T k H J \k B k Q k H J k h H 

Here, H and h are number densities associated various combinations of the up- and down- 
type Higgs supermultiplets in the MSSM (defined below); is the corresponding vector 
current density; Q and (B,T) are the number densities of particles in the third genera- 
tion left- and right-handed quark supermultiplets, respectively; the kn,h,Q,T,B are statistical 

weights; S^f is a CP-violating source; and T H , T Y , f y, and T Y are transport coefficients. 

Physically, the presence of results from an imbalance between the rates for particle 
and antiparticle scattering off the bubble wall, favoring the generation of non- vanishing su- 
permultiplet densities H and h. In contrast, the terms proportional to T# and Py cause 
these densities to relax to zero. The terms containing Ty and Ty favor chemical equilib- 
rium between Higgs superfield densities and those associated with quark supermultiplets. 
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To the extent that the rates Ty and Ty are fast compared to the rate of relaxation, any non- 
vanishing Higgs supermultiplet density quickly induces non-vanishing densities for quark 
supermultiplets, thereby facilitating EWB. Understanding the microscopic dynamics of this 
competition between CP-violating sources, relaxation terms, and Higgs-baryon sector cou- 
plings is essential to achieving a quantitative description of EWB. 

In previous work, we computed Th and using the Closed Time Path Schwinger- 
Dyson equations and considering the lowest-order couplings between superfields and the 
spacetime varying Higgs vevs. Here, we focus on the terms proportional to Ty, Ty, Ty that 
are generated by Hqq Yukawa couplings, the corresponding supersymmetric interactions, 
and the SUSY-breaking triscalar couplings 1 . We make several observations regarding these 
terms: 

(i) In previous treatments, Ty and Ty were estimated from scattering processes such as 
tn + g — > ti + H®, making them proportional to one power of the strong coupling, a s . 
We find, however, that there exist contributions to Ty occurring at zeroth order in a s 
that are generated by decay and inverse decay processes such as £r + £l <-> H®. To the 
extent that the three-body processes are kinematically allowed, their contribution to 
Ty can be considerably larger than those generated by scattering. We also show that 
|ry/IY| is typically < 1/2 for MSSM parameters consistent with precision electroweak 
data and the existence of a strong first order phase transition. (The authors of Ref. [Io| 
argued that |fy| = |Ty|.) We solve the transport equations numerically and find that 
inclusion of the three-body contributions affects the baryon-to-entropy ratio Yb at the 
10 — 20% level for realistic choices of the MSSM parameters. We provide a detailed 
analysis of the dependence of Yb on Ty and the MSSM parameters that determine it. 

(ii) In most of the early studies of EWB in the MSSM, it was assumed that the rate 
Ty of Yukawa-induced processes is "fast" compared to all other relevant time-scales, 
implying that the Yukawa-induced transfer of non-zero Higgs/Higgsino density to non- 
vanishing chiral charge density is more efficient than relaxation. This assumption has 
motivated an expansion in powers of 1/Ty. We show that there exist corrections to 
the Higgs density at linear order in this expansion that have not been included in 
previous treatments. After including these terms, we find that the expansion itself 
breaks down - even for the enhanced values of Ty that result from inclusion of the 
three-body contributions - due to the presence of chirality-changing processes in the 
bubble wall whose rates Th and Tm can be larger than Ty. We study numerically the 
impact of keeping a finite Ty: we find that the corrections to the Ty — > oo limit of Yb 
range between 20% and 100%, depending on the values of the other rates. 

(iii) The terms containing Ty and f y have been included in the earlier studies of Refs. 0, 
[ill I20L EH ] 2 , whereas the one involving fy is new. In the MSSM, one often assumes 
that the triscalar coupling involving the down-type Higgs scalars, the doublet scalars 
Q, and the right handed scalars b is proportional to the bottom Yukawa coupling, ?/&. 
For tan/5 ~ 0(1), one has Ub/yt <C 1 and the impact of the Ty term is relatively 
minor. For scenarios with large tan/3, however, need not be small compared to 



1 In previous work, only the contributions to terms of this type generated by Standard Model Yukawa 
interactions were considered, leading to the use of the subscript ll Y" . 

2 In the notation of Ref. Q, fy = -pT Y . 
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y t . In this case the transport coefficient Ty and other terms (not shown) that couple 
to the B supermultiplet need not be suppressed, and the coupled set of transport 
equations must be augmented to include dynamical 6-quarks and their superpartners. 
Although in the present study we do not consider this large tan (3 scenario, we provide 
the general formulas that allow one compute Ty. 

In the remainder of the paper, we discuss our detailed analysis of the Ty-type terms 
that lead to these observations. In Section |H1 we consider these terms for generic Yukawa 
and tri-scalar interactions and analyze their dependence on the relevant mass parameters. 
In Section HID we specify to the MSSM, including detailed analytic and numerical studies. 
Here, we include contributions from both SM particles and their superpartners (in contrast to 
previous anlayses that included only SM scattering terms), and note that the superpartner 
contributions tend to increase the magnitude of Ty. In Section IIVI we solve the coupled 
transport equations to obtain the baryon-to-entropy ratio, and show why one would not 
expect an expansion in 1/Ty to yield a reasonable approximation to the exact solution. We 
summarize this work in Section |V| Various technical points are discussed in the Appendices. 

II. THREE-BODY SOURCE TERMS: BUILDING BLOCKS 

Our approach for deriving the source terms in the quantum transport equations is based 
on the Closed Time Path Schwinger-Dyson equations. An extensive discussion of this frame- 
work is given in our earlier work [lij . Here, we give a brief summary of our method and 
use it to derive the source terms generated by supersymmetric Yukawa and SUSY-breaking 
tri-scalar interactions to leading order in the loop expansion. 

A. Formalism and method 

Ordinary quantum field theory is not appropriate for treating the microscopic dynamics of 
the electroweak phase transition (EWPT), since the non-adiabatic evolution of states and the 
presence of degeneracies in the spectrum break the zero-temperature, equilibrium relation 
between the in- and out-states. The non-adiabaticity arises because particle interactions 
occur against a spacetime- varying background field (the Higgs vevs), while thermal effects 
associated with non-zero temperature introduce degeneracies in the spectrum. The impact 
of non-adiabaticity and degeneracies on quantum evolution can be treated systematically 
using the Closed Time Path (CTP) formalism j^]. In this formulation the time arguments 
of all fields and composite operators lie on a path V that consists of a positive branch V + 
from — oo to +00 and a negative branch P_ running back from +00 to —00. Fields whose 
arguments lie on V+ precede those on V- along the path V . Moreover, those lying on V + 
are time-ordered while those on V- are anti-time-ordered. 

With this prescription the standard time-ordering operator T is replaced by the path- 
ordering operator T-p and the perturbative expansion is formally identical to the equilibrium 
case. In applying Wick's theorem, however, one must allow for contractions involving all 
possible combinations of fields taken from either V+ and V-, leading to a generalized Green's 
function that accounts for path-ordering. Specifically, the bosonic and fermionic Green 
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functions are given by 



G(x,y) = (T P [0(x)0%)]> 
S(x,y) = (T v fy(x)li{y)]) 



(2) 
(3) 



where (...) denotes an average over the physical state of the system, which may be described 
by an appropriate density matrix. In practical applications it is convenient to use ordinary 
time arguments, in terms of which each of Eq. (j2J) and (j3J) represents four Green functions 
and decomposes in various components. To establish the notation we recall here explicitly 
the bosonic Green functions: 

(4) 
(5) 
(6) 
(?) 

where the superscripts "a" and "6" in G ab (x,y) indicate the branch V± on which the time 
components of x and y lie, respectively and where T is the anti-time-ordering operator. 

The equations governing the spacetime dependence of number densities of a given bosonic 
or fermionic species can be derived from the Schwinger-Dyson equations for the generalized 
Green's functions G(x,y) and S(x,y) and have the following form |qL J23t] : 
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E>(X, z)S < {z, X) - S>{X, z)Z<{z, X) 



+S<(X, z)T>{z, X) - S<(X, z)S>(z, X) 



(9) 



The RHS involves a causal time integral over the system's history and is expressed in terms 
of the Green functions (J2J), (JHJ) and self energies T*b,f that encode all the information about 
particle interactions. This feature allows for a consistent treatment of both CP-violating 
terms "sourcing" a given particle density as well as CP-conserving interactions that tend to 
transfer this density to other species or cause it to relax away. Previously, the leading CP- 
violating contributions t o E n F gen erated by scattering from the Higgs vevs (see Fig. |TJ) were 
computed in Refs. 0, lllUl5lTl7l| while the corresponding CP-conserving relaxation terms 



generated by the same processes were derived in Ref. [15j- Here, we extend these analyses 
to include the three-body source terms that arise from Yukawa and tri-scalar interactions 
at one- loop order (see Fig. |2J). 

In general, the Green's functions (J2J), © are dynamical objects that can be obtained 
by solving the transport equations (jSJ) and Q- However, the hierarchy of time and energy 
scales present during the electroweak phase transition allow for simplifications in treating 
the transport equations [l5j]. The time scales are a decoherence time, Td, associated with the 
departure from adiabatic evolution; a "plasma" time, r p , associated with mixing between 
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FIG. 1: Leading contributions to the self-energies T,b,f generated by scattering from Higgs vevs. 

degenerate states in the finite temperature spectrum; and the "intrinsic" quasiparticle evo- 
lution time, Tint, associated with time evolution of a state of definite energy. In terms of 
physical parameters associated with the plasma, one has Td ~ l/(v w k e g), r p ~ 1/T P , and 
Tint ~ 1/^5 where v w is the velocity of expansion of the bubble wall; k e s is an effective 
wave number that in general depends on the quasiparticle wave number and wall thickness, 
L w ; T p is the thermal width of the quasiparticle; and uj is the quasiparticle frequency that 
depends on both the particle momentum and thermal mass. For the EWPT, one has that 
£d = T int/ r & <C 1; £ p = T int /r p <C 1; and e^/sp = T p /r d <C 1. In addition, the small densities 
present at the EWPT imply a hierarchy of energy scales: ^ = /i/T C 1, where \x refers to 
the chemical potential of any particle species. The existence of these hierarchies allows for 
a number of simplifying approximations in solving the transport equations: 

(i) Use of the quasiparticle ansatz for the Gi(x, y). This relies on e p < 1, that is, that the 
damping rates T p that broaden the spectrum of excitations are typically suppressed 
when compared with the excitation frequencies {Ti/uji <C 1). 

(ii) Working near kinetic and chemical equilibrium. This approximation relies on T p /Vd <C 
1, that is, the plasma interactions among quasiparticles are fast compared to the 
decoherence time, thereby leading to approximate, local equilibrium among quasipar- 
ticle species. Consequently, one may approximate quasiparticle distribution functions 
appearing in the Green functions by their equilibrium forms and track quasiparticle 
densities with local chemical potentials. The error engendered by doing so is 0(e^/e p ) 
and is, thus, negligible. 

Motivated by these considerations we evaluate the source terms on the RHS of the trans- 
port equations (JHJ) and (jHJ) using the free-particle form of the Green functions. For example, 
for the boson Green functions, we have 

Pl (k 0l k) (10) 

Gf(x,y) = J -^-e-^-^fBiko,^) Pi (k ,k) (11) 

with spectral functions Pi(ko, k) = it/u^ [5(k° — cj k ) — 5(k° + (cj k = Vk 2 + m?) that 
can be appropriately modified to take into account collision-broadening and thermal masses, 
and distribution functions close to the equilibrium form 

fB{ko,Hi) =n B (k ,iii) +0(e d /e p ) , (12) 




UK 

(2nY 



-ik-(x—y) 



1 + /b(&o, lh 
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where ^(fco, aO = l/^ k °~^ i ^ T — 1] and ^ is a local chemical potential . 

Upon expanding the source terms to lowest non-trivial order in the £d,p,^ and relating 
current and chemical potential to local densities through 

• h {x) = -A Vm(x) m{X) = ^ Hm.lT) ^(x), (13) 

o 



where k^rrii/T) is a statistical factor (see, e.g. |15j). we obtain the quantum transport 



equations 3 



h i -D i V i n i = S i [{n j }] . (14) 

In Eq. (|14|) both CP-violating effects and relaxation rates are encoded in the quantum 
mechanical sources ^[{riy}]. 



B. Results for generic tri-scalar and supersymmetric Yukawa interactions 

Let us consider now the generic three-scalar interaction, 

Ant = KA4>l<P*r4>h + h.c. , (15) 

where X s is a dimensionless coupling and A s is a mass scale 4 . This interaction generates 
contributions to the self-energy appearing on the RHS of Eq. (JHJ) through the one-loop 
diagram depicted in Fig.^c). As an example, we give the self-energy for the complex scalar 

^' j/) = -\XsA s \ 2 G > L ' < (x,y)G > H ' < (x,y) , (16) 

Importantly, the RHS of Eq. (|TB|) is manifestly independent of possible CP-violating phases 
appearing in the coupling X S A S and therefore does not contribute to the CP-violating source. 
We obtain similar results for the self-energies of (pt and ipn- This situation contrasts with 
that for the Higgs vev scattering contributions derived from Fig. ^ where interference terms 
involving the up- and down-type Higgs vevs at the different vertices contain CP-violating 
phase effects. 

Inserting Eq. (|16|) into Eq. (jSJ), using the Green Functions of Eqs. (jl(Jlllll2j) . and ex- 
panding to first order in and zeroth order in e^ p (setting the thermal widths to zero), we 
obtain the leading, three-body contribution to Sr on the RHS of Eq. (|14j) . We find that the 
three-body sources for the particle number densities of the complex scalars (ft l, (fin and 
are related to each other and are given by 

S R {X) = -S L (X) = -S H (X) 

= - (Vr - PL - ^h) (X) \\ s \ 2 1 B (A S ; m R ,m L ,m H ) , 

3 The quantum transport equations (|14|) are sometimes referred to as quantum Boltzmann or diffusion 
equations. 

4 In the MSSM, A s is the Yukawa coupling and A s is either the /x-parameter or the soft, tri-scalar coupling. 



8 



...ry. 



FIG. 2: Self energies for scalar and fermion fields induced by the Yukawa and triscalar interaction 
lagrangians of Eqs. (|15|) and 1|2U|) . 



in terms of the function 



T B {A S ; m R , m L , m H ) = - ^^j, J duj R j 



\Af f°° f u i 



x jr^u^) [1 + n B (uJ L )]n B (uJ L - u R ) [6{m R - m L - m H ) - 9(m L - m R - m H )] (18) 
- n B {uj R )n B {uj L ) [l + n B {u L + Lo R )]d(m H - m R - m L )j, 



with integration limits given by: 



1 

2m 



,2 



u l = KTT^R \ m R + m L~ ™ H \ 



± [(u% - m 2 R ) (m 2 R - (m L + m H f) (m 2 R - (m L - m H ) 2 )] 



1/2 



(19) 



The presence of mass thresholds and combinations of Bose distributions in Eq. (|18|) 
makes clear its interpretation in terms of physical processes in the plasma: decay R — > 
L + H and all possible emission/absorption channels that are kinematically allowed. It 
is straightforward to integrate over uol and obtain a representation of the source in terms 
of one-dimensional integrals. We give this formula in Appendix [X] Finally, we note that 
Sr,l,h are of first order in the e counting discussed above, whereas the leading CP-violating 
sources and CP-conserving relaxation terms generated from the tree-level graphs of Fig. ^ 
are 0{e^e v ) and 0(e p e^), respectively. Nonetheless, they can be similar in magnitude to 
S r ,l,h since the latter contain additional phase-space suppression factors ~ 16n associated 
with the absorptive part of one-loop graphs. 

We now consider contributions from a generic Yukawa interaction 

Ant = \f{UlP L ^2 + (20) 

that generates contributions to both scalar and fermionic self-energies on the RHS of Eqs. (jSJ) 
and through the diagrams depicted in Fig. I21^a),(b). The resulting source for particle 
number densities associated with the complex scalar and Dirac fermions ipi and fa are 
related to each other and read: 

S<j>{X) — -S^^X) — S^ 2 (X) 

= - (/^ - A*i + // 2 ) (X) \X f \ 2 l F (m 1 ,m 2 ,m (j) ) , 
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FIG. 3: Left panel: Tp/T 3 as a function of ms/T for mi/T = 1 and m 2 /T = 1, 2, 5. Right panel: 
Ib/T 3 as a function oimu/T for j4 s /T = 1, uil/T = 1 and mji/T = 1,2,5. 



where 

l F (m 1 ,m 2 ,m 4 



16tt 3 T 



(m{ + 



//?., 



gL>i 



x < riBiu)^) [l - n f (wi)]n F (wi - u^) [9^ - m 2 - m$) - 6{m^ - m 1 - m 2 )] ( 22 ) 



+ [l - n F {u x + ^)]6>(m 2 - m x - m^) 

with integration limits on cj^ given by: 
1 



2m( 



ui + m l 



— m n 



± [(ul - ml) (ml - (m 2 + m^) 2 ) (m\ - (m 2 - m<p) 2 )] 



1/2 



(23) 



As in the bosonic case, Eq. (|22|) has a direct interpretation in terms of decay, emission and 
absorption of 0, ipi and ip 2 in the plasma. The integration over is straightforward and 
we report its result in Appendix 1X1 

Eqs. (I17H18JI and (j21H22J) are central new results of this paper and represent the building 
blocks out of which we can construct the three-body physical sources in the MSSM (Section 
TTTj) . In order to identify the dominant contributions to the MSSM sources, where many 
individual building blocks contribute, it is instructive to characterize the behavior of 1b,f 
as function of the masses of the interacting particles. The main features are: 

i) Symmetry properties under exchanges m^ <-» m R ; «-> m B and m\ <-> m 2 : 

Z B {A S ; m R , m L , m H ) = X B {A S ; m L , m R , m H ) If {mi, m 2 , m<j) = X F {m 2 , m\, m<$) 

ii) There are threshold effects which can be read off via the explicit ^-functions. In order 
for the rate to be non-zero, the mass arguments have to be such that at least one of 
the two body decays a — > b + c is kinematically allowed. 
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iii) Ib{A s ; m R , ttil, tuh) and Xp(mi,m2,m^,) are largest when the three masses are such 
that the largest mass is slightly greater than the sum of the two smaller ones (just 
above threshold). 

iv) Tb(A s ; mn,mL,mH) and If(ttii, m,2, m^) become vanishingly small as any of the 
masses becomes much larger than the temperature. This reflects Boltzmann sup- 
pression of the thermally averaged rate. Moreover, Xp(mi, rri2, m^) vanishes as all the 
masses become much smaller than the temperature. 

The above properties are illustrated in Fig. |31 where we plot Ip/T 3 as a function of m^/T 
for representative choices of mi and m 2 (left panel) and similarly X R /T 3 as a function of 
vTLujT for representative choices of m^, m^, and A s (right panel). 

III. THREE-BODY SOURCE TERMS IN THE MSSM 

The results of Sec. HU allow us to calculate the sources for quark, squark, Higgs, and 
Higgsino particle densities generated by the supersymmetric Yukawa and SUSY-breaking 
triscalar interactions in the MSSM. We focus on those involving the third-generation quark 
supermultiplets whose interactions generally depend on the large Yukawa coupling y t . As 
noted in the Introduction, these interactions dominate for tan j3 ~ 0(1), whereas interactions 
proportional to ?/& can be important for large tan j3. While the results in the previous section 
would allow us to compute these y^ effects - such as the transport coefficient Ty appearing in 
Eq. (JTJ) - including them would lead to a more complex set of coupled transport equations. 
For simplicity, we focus here on the smaller tan (3 case with a/ oc yf - wherein interactions 
involving y t dominate - and defer a more general treatment to a future study. 

A. Interactions in the MSSM 

The terms in the MSSM superpotential generating interactions proportional to y t are: 



where the weak doublets are defined Q3 = (ti, H u = H®), and Hd = (H d , H d ) In 
addition the soft SUSY-breaking Lagrangian contains the terms: 



In the minimal supergravity (mSUGRA) scenario for SUSY breaking, the a-parameters are 
proportional to the Yukawa couplings, e.g. a t = ytA t for some mass parameter A t . Thus this 
part of £ so f t also generates contributions to the top three-body source that are proportional 
to y t . 

From both the supersymmetric and soft SUSY-breaking sectors, we obtain the tri-scalar 
interactions 



W = y t Q 3 H u t R + fiH u H d , 



(24) 



£ so ft = —(itQ?,H u t* R + h.c. 



(25) 



Y 

scalar 



y t t* R t L (A t H° u + fi*H° d *) + y t t* R b L (A t H+ - //#"*) + h.c, 



(26) 



and the supersymmetric Yukawa interactions 



Y 

fcrmion 



y t (-H%t L + H+t R b L ) + y t (-t R t[H°J + i R b{H+J) 
+ y t (-t L t R H° u + b L t R H+) + h.c. 



(27) 
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In order to write this Lagrangian in the form appearing in Eq. (|2U|). we combine the two- 
component Higgsino spinors into four- component Dirac spinors, which is sensible in the 
unbroken electroweak phase where the mass terms for Higgsinos are simply: 



£ mass = -fiH+Hj + pH$H° + h.c. (28) 

First rotating the fields H® ,+ — > e - "^ H®' + to remove the complex phase from fi, we define 
the Dirac spinors 

H + \ ~r f-H G . 



which have Dirac mass We define chemical potentials / u, i ^+,/i^- corresponding to the 
vector charge densities \j>7°\l/ for these Dirac fields. In terms of these fields, the Yukawa 
interaction terms are: 

£Lnio„ =yt(-H°j R p L t L + H+i R p L b L ) 

+ yte^{l R i L PR^% + t R b L P R ^% + ) (30) 
+ yte-^itLtuPL^HO + b L t R PL^H+) + h.c, 
making use also of the charge-conjugated fields: 



= (lit) *S» = (St ) . 



where ^ c = C^ T , with C = ry 2 7° 



B. Source Terms in the MSSM 



Having identified the relevant interactions in the MSSM Lagrangian proportional to y t in 
Eqs. (|2*Bj) and (j3*Uj) . we can write the sources for the densities of the particles appearing in 
these interactions using the general results of Eqs. ()17I21|) . 

To be concrete, let us focus on the right-handed top squark and quark densities. Similar 
formulas will hold for the left-handed squarks and quarks, and the Higgs and Higgsinos. 
The source for the right-handed top squark number density rif is: 



S U X ) 



N cVt [(lJ-i R ~ M h ~ ^Ho)I B (At] m iR , m h ,m H o) 
+(m r ~ M L ~ ^H+) x B{A t - m iR , m~ bL ,m H +] 

+(m r - m l + hh*Pb{h\ m R ^i h ^Hi) 

+ (m r - M L + /%- ) x b {h; m~t R , m~ H , m H - ) 
+(M R - Ml ~ ^Ho) X Firn n+} m tL ,m iR ) 
+(M R ~ Ml ~ VH+) X F('rri il0} m bL ,mi R ) 



(32) 
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and for the quark density n tR : 



R 

S t R ( X ) = -N c Vt [(lh R ~ IH L ~ ^Ho)T F (m tR ,m tL ,m H o) 
+ (Vt R ~Vb L ~ fi H +)l F {m tR ,m bL ,m H +) 
+ (Vt R ~ Pt L ~ VHo)lF(m tR ,m m , m ih ) 



mi 



The various chemical potentials appearing in the source can be related by making the 
assumption, first introduced in Ref. |4j, of fast gauge and gaugino interactions and zero 
density of gauge bosons or gauginos (/iy = fly — 0). In this case, pairs of superpartner 
densities are in chemical equilibrium, as are members of the same gauge multiplet. Thus, 

th B = m R = (34a) 

Vt L = m L = Vb L = H L = ( 34b ) 

= /%+ = Hh u (34c) 

^ H ° d = = m A (34d) 

Hh+ = ^ho = t- L H (34e) 

Relating the scalar Higgs chemical potentials fJ*H ud to the Higgsino chemical potential [i^ is 
somewhat more subtle and we refer to Appendix Q5] for a derivation. Defining the combina- 
tions, 

fiH = g (Ph v ~ VH d ) (35) 
l*h = \ (HH u + VH d ) , (36) 
the supergauge equilibrium condition reads: 

t*H = ■ (37) 

As noted in previous work, the assumption of supergauge equilibrium - together with the 
relations I3*T|) - suggest combining the various particle densities in equilibrium with one 
another into: 

T = n tR + n- tR (38a) 
Q = n tL + n bL + n h + n~ bL (38b) 
H = n H + + nijo — n H - — n H o + rif I+ + ngo (38c) 
h = n H+ + n H° + n H - + n H o (38d) 

Adding together the top and stop sources in Eq. (J32I33|) and using the relations among 
chemical potentials 1331 l3~7j) leads to the Yukawa source for the density T reported in 
Eq. (jBTl|) of Appendix [Bj Finally, by noting that the masses of weak doublet partners are 
the same (see Eq. (|B12jl ) and converting the chemical potentials to densities using Eq. (fT3*j) 
we obtain 

S*(X) = -Ty(L -f)-r Y ± (39) 

\k T k Q k H J k h 
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where 



r 



Y 



J-2 



r 



y 



T B (A t ; m~ tR i nig, m Hu ) + J B (/x; m tfl , m^, m^J 
+J F (/i, mg, m t -J + l F (m tR , m Q , m Hu ) + l F (m tR , tUq) 

12 ^2 Vt Xb ^ m t R i m & m H d ) ~ ?b{Au m iR ,rriQ, m Hu ) - Z F {m tR ,m Q , m Hu ) (40b) 



Similar formulas hold for the sources Sq h h . 

C. Transport Equations and study of Yukawa rates 

Incorporating the Yukawa contributions to the sources into the full set of transport equa- 
tions for the densities T, Q, H derived in Ref. [l5j . we obtain 



"~ M \k T knl M \k T knl 1 



TQH\_ T h +v j2QT^ { Q + T) 



k T k Q k H J k h \k Q k T k B 
Q T\ .(T Q 



(41a) 



ry[ Q_L + IL) + f Y ±-2T ( 2Q T ■ 9 ^ + T > 



kQ k T k H J k h \k Q k T k B 



(41b) 



PH, = -r E ^-r Y (£-^ + P\ + r Y ± + s¥ (4ic) 



M k H \k Q k T k H ) k h H 

In addition, there should be one more equation for d^h^, but we have left for future work 

cyp 

the calculation of the CP- violating contribution, S{ to its source, as well as its relaxation 

h 

coefficient Y h . 

The structure of the transport equations (}4~Tj) is similar to that of the equations derived 



transpi 



in the treatment of Ref. 14, llQl 11 . However, use of the CTP framework leads to a number 



of new features that we highlight: 

i) The appearance of new combinations of densities that do not arise in earlier treat- 
ments - such as those involving - follows from a systematic treatment of the CTP 
Schwinger-Dyson equations. 

ii) The Yukawa rates Ty and T Y arise at lower-order in a s than the corresponding terms 
in previous treatments. As indicated in the Introduction, these rates were calculated 
to 0(a s ) from scattering processes such as t R + g — > t L + H® and only the contribu- 
tions from Standard Model particles were included. We have included here the O(a s ) 
contributions generated by decays and inverse decays within the plasma, which - when 
not vanishing due to threshold effects - can be of comparable size or larger than the 
0(a s ) scattering terms. This can be appreciated by comparing the behavior of Ty 
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from decays (Eq. PU|)) and from scattering (see Ref. 



r 



Y 



i scattering 



idecays 
Y 




8T 2 



2 q (T) 



) 



6ir 3 



VtT x 0(1) , 



(43) 



(42) 



where M is a typical (thermal) mass of the order of the electroweak scale (could be a 
soft SUSY breaking mass term), m q (T) is the quark thermal mass, and (3 = 1.202. 

iii) Because we have included both SM particle and superpartner contributions, Ty and Ty 
display a non-trivial dependence on the MSSM par ameters. Similar observations have 
been made about the CP- violating sources P. lid llif and leading chiral relaxation 
terms [l^|, for which the possibility of resonant enhancements have been observed. 
We note that the enhancements of the CP-violating sources and chiral relaxation 
are generally not accompanied by resonant enhancements of the Ty and Ty terms, 
thereby leading to a more subtle competition between the effects of CP-violation, 
chiral relaxation, and density transfer. 

A quantitative illustration of the above points ii) and iii) is given in Fig. HJ where we plot 
Ty and Ty versus the MSSM parameter \ u\ f or T = 100 GeV. In the numerical evaluation we 
include thermal masses as calculated in [24j and, for illustrative purposes, we use the weak- 
scale SUSY parameters given in Table |l] consistent with electroweak symmetry breaking, a 
strongly first-order electroweak phase transition and electroweak precision tests. The non- 
trivial \i dependence displayed by Ty cay is due to threshold effects in the functions T B F . 
The dashed straight line in Fig. 0] represents r^ atterm s j n l ar g e regions of parameter space 



We conclude this Section by noting that for typical values of SUSY parameters, the chiral 
relaxation rates Tm,h (active only in the broken electroweak phase) are of comparable size 
or larger than Ty. All of these rates, in turn, are much larger than the diffusion rates v^/Di, 
which for typical values of the diffusion constants 0, 0] and wall velocity vary in the range 
10 -3 — 10 -2 GeV. We discuss the consequence of this when solving the diffusion equations 
in the following Section. 

IV. SOLVING THE TRANSPORT EQUATIONS AND PHENOMENOLOGY 

The baryon asymmetry is seeded by the density of left-handed weak isodoublets ul = 
5Q + 4T p, 0], which we obtain by solving the transport equations (JHJ). In this section 
we study the impact of Ty on the solution of the system (|41j) and on the overall baryon-to- 
entropy ratio Y B = n B /s. 

Before entering the details of our analysis, let us shortly recall the basic notation (see 
[l5| and references therein) and describe the input MSSM parameters which will be used in 
the subsequent numerical explorations. The baryon-to-entropy ratio can be expressed as an 
integral of ul = 5Q + 4T in the unbroken phase: 



j2 1 T-\dccay . t-vi 

we find 1 v > 1 



1 scattering 




(44) 
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100 GeV and SUSY mass parameters as described 



in the text. In large regions of parameter space we find T Y 



decay 



-pscattcring 
> ly 



where T ws is the weak sphaleron rate T ws = Qncc^T (with k ~ 20 |25j), rip is the number of 
fermion families, D q is the quark diffusion constant, v w is the wall velocity and 



n 



(v w ± \ f' 

2D„ \ V 



vl + 4D a K 



^squark \ ^ 

6 / 



(45) 



where n squar k is the number of flavors of light squarks. Isolating the dependence on the 
CP- violating phases 0„ and 4>a, Yb is conveniently parameterized as follows (T - " 



Y B = F 1 sin 0„ + F 2 sin i 



(46) 



in terms of Fi (arising from the Higgsino source) and F 2 (arising from the squark source). 

In all the plots reported in this section, we adopt for the weak-scale SUSY parameters 
the reference values reported in Table HJ which are consistent with a strongly first-order 
electroweak phase transition and the constraints from precision electroweak physics as well 
as direct searches. Note that a CP-odd Higgs mass ttla = 150 GeV translates into A/3 ~ 
0.015 |26|| . From the reference values of Table H] one can derive typical values for the bubble 
wall velocity and thickness, for which we use v w = 0.05 j2?| and L w = 25/T j26{. With 
this choice of parameters one has F 2 ~ lO -3 ^. 

We now discuss in greater detail the role of Yukawa-induced rates on the transport equa- 
tions. 



A. Revisiting the approximation of fast Ty: need for numerical solution 

Starting with the work p| , the conventional practice has been to solve the system of trans- 
port equations (|4T|) under the assumption that the rate T Y of Yukawa-induced processes 5 



As well as the rate T ss of strong sphaleron processes. 
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tan (3 = 10 
M iR = 
Mq = 1 TeV 
M 2 = A t = 200 GeV 
m 2 H = -(100 GeV) 2 




100 GeV <\fi\< 400 GeV 



TABLE I: Reference values of weak-scale SUSY parameters. 



is "fast" compared to all other relevant time-scales, thereby ensuring a chemical equilibrium 
condition among H, Q, and T. Doing so allows one to obtain analytic expressions for Yg. 
The assumption of fast Yukawa interactions is well justified in the unbroken phase ahead of 
the advancing bubble wall, where a particle may diffuse for a period characterized by the 
the inverse of the diffusion rate T^is = v 2 w jD before the bubble wall catches it. In order for 
Yukawa processes to be effective in this region, they must act quickly on the time scale r^L, 
and one, indeed, finds that Yy ^> Y^is for typical values of the diffusion constants |H 
and wall velocity. In the broken electroweak phase, however, Y^ s is no longer the only 
relevant time scale. In addition, Yukawa processes must compete with scattering from the 
spacetime- varying Higgs vevs that leads to relaxation of chiral charge and Higgs supermulti- 
plet densities. Importantly, the corresponding rates (Y M and Y H , respectively) are as large 
as or larger than Yy - even after including the O{a s ) contributions to Yy. As a result, the 
interplay of these competing processes within the bubble wall is significant, and imposing 
the condition of Ty-induced chemical equilibrium is not justified. 6 

To make this key point more explicit, we have solved the transport equations in powers 
of l/ry jSS and analyzed the magnitude of the corrections to the Ty^s — * oo limit. Explicit 
details are given in Appendix O where we point out that the most important correction was 
missed in previous analytic approaches to this problem - namely the correction to the H 

density induced by an effective shift in the source S^. The analysis of Appendix [0 implies 
that fractional corrections to the baryon asymmetry to first order in 1 /Yy read 



where rp ~ 0.07 . Substituting the earlier estimates of Yh and Yj d |4j into this expression, 
we find 5Yb/Yb ~ 0.1 - indeed a small correction. However, when using Y^ and Yy 
as calculated in Ref. and the present work within the CTP framework, we find much 
larger corrections: 5Yb/Yb ~ This difference is due primarily to the larger values 

of IV and Ym obtained in our framework 0] (even off resonance) compared to previous 
calculations 0, H, lf|. 



6 In addition, the authors of Ref. |20j noted that the condition Ty — * oo causes a parametric suppression 
of the Higgs source h, while for realistic parameter choices, the suppression factor turns out to 0(1). 




(47) 
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FIG. 5: We plot here the ratio F 1 /Yg mAP versus Ty for two values of the SUSY parameter: 
|/i| = 200 GeV (solid line) and = 250 GeV (dashed line), corresponding to on-resonance and 
off-resonance baryogenesis, respectively. All other parameters are fixed at the reference values of 
Table E We use the central value yWMAP = Q 2 x 1Q -n Q 




The above considerations imply that in order to avoid 0(1) uncertainties in the calculation 
of Yb-, one requires a full numerical solution of the system (|41|) . In order to quantify the 
effect, we plot in Fig. El the ratio Fi/Yg fMAF versus Ty for two values of the SUSY fi 
parameter: \f/,\ = 200 GeV (solid line) and = 250 GeV (dashed line), corresponding to 
on-resonance and off-resonance baryogenesis, respectively. All other parameters are fixed as 
in Table HJ Typical values of Ty lie in the range 5 — 10 GeV (see Fig. HJ). The curves in 
Fig. El illustrate two key points of the Yukawa- induced dynamics: 

(i) Efficient chargino/neutralino- mediated baryogenesis occurs for Ty 3> v^/D^ ~ 0.0025 
GeV, as the Higgs supermultiplet density H injected in the unbroken phase is efficiently 
converted into LH top-quark density (fueling sphaleron processes) before the bubble 
catches up. Inclusion of the 0{pP s ) terms in Ty affects Yb at the 10 — 20% level, as 
one is already in the plateau region in Fig. El 

(ii) As Ty increases (keeping all other rates fixed) the baryon asymmetry reaches a max- 
imum and then starts decreasing towards its asymptotic value. This behavior can 
be understood qualitatively as follows. In the non-resonant case (dashed line), as Ty 
increases, Yukawa induced processes start to complete with Th inside the bubble wall, 
thereby transferring H density to Q, T densities. The latter subsequently relax away 
due to processes or diffuse very inefficiently into the unbroken phase. This effect is 
less pronounced in the resonant case (solid line), where Yb first grows as Ty becomes 
more efficient compared to diffusion ahead of the bubble wall, but then saturates due 
to the presence of resonantly-enhanced Higgs supermultiplet relaxation within the 
plasma. 

Summarizing, the main message emerging from Fig. El is the following: keeping Ty finite 
and in the realistic range (few GeV) can increase Yb by a factor between 20% (resonant 
case) and 100% (nonresonant case) compared to the Ty — > oo limit. 



18 



100 




100 150 200 250 300 350 400 
11 (GeV) 

FIG. 6: F\ (solid line) and F% (dashed line) versus with all other parameters fixed at the 
reference values of Table [TJ 

B. Phenomenology update 

In a consistent analysis Ty should not be treated as an independent quantity (as we did in 
the last section for illustrative purposes) but rather as a function of the MSSM parameters 
(as we did in Section ITTT|) . Doing so after numerically solving the transport equations, we 
study the behavior of (Eq. (jH2)) as a function of the MSSM parameters. For illustration, 
we show in Fig.|H]the dependence of F\ (solid line) and F2 (dashed line) on with all other 
input as in Table UJ The plot highlights the resonant behavior of F\ discussed in 0, 0] . 
The behavior of F2 follows from the fact that F2 is proportional to |ju|(r# + T^) -1 / 2 : the 
dip at ~ M% ~ 200 GeV reflects the resonant enhancement of IV. The overall scale of 
Fx 2 is set by A/3 which in turn depends crucially on the CP-odd Higgs mass tjia [13 : here 
we use rriA = 150 GeV but one should keep in mind that higher values of can lead to 
sizable suppression of 

Finally, we investigate the impact of electric dipole moment (EDM) searches on this 
particular EWB scenario. It has long been recognized that, given the spectrum of super- 
symmetric particles, constraints from the electron j2^|, neutron |3(|, and nuclear (3j| EDMs 
pose tight limits on the size of CP- violating phases (for a review see 0). These could ul- 
timately enter in conflict with the requirement of successful baryon asymmetry generation, 
making EDM searches a great discriminating tool among theories of baryogenesis. 

To illustrate this point we plot in Fig. [7| the allowed bands in the 0^-0^ plane resulting 
from present limits on electron and neutron EDMs and successful baryogenesis, for a given 
choice of the SUSY mass parameters. Here, we have employed one-loop SUSY contribu- 
tions |32|]. We take the first- and second-generation sfermions, as well as the gluinos, all 
degenerate at 1 TeV, while all other input is fixed as in Table |1J In the left-hand panel 
we use M 2 = \fj,\ — 200 GeV (resonance peak), while in the right we use M 2 = 200 GeV, 
\n\ = 250 GeV. 

Figure [7| illustrates the complementarity of various EDM measurements in constraining 
the new CP-violating phases in general. It also shows that in this particular scenario it is 
the electron EDM that poses the strongest constraints on electroweak baryogenesis. In order 
to quantify the dependence of the EDM constraints on the heavy sfermion masses, we plot 
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FIG. 7: Allowed bands in the fypr^A plane implied by consistency with the 95% C.L. limits on 
electron and neutron EDMs and baryogenesis. The solid lines correspond to the constraint from 
the electron EDM (|cZ e | < 1.9 x 10 _27 e • cm [ioT]). and the dashed lines correspond to the neutron 
EDM (|cZ n | < 3.6 x 10~ 26 e • cm HU). These EDM constraints correspond to sfermion masses (mj) 
fixed at 1 TeV. The shaded EWB band is the reg ion consistent with Yb from BBN jsif at 95% C.L. 
(which includes the Yb range from WMAP 28]). In the left-hand panel we use \fj,\ = M2 = 200 
GeV (resonance peak), while in the right-hand panel we use M2 = 200 GeV and = 250 GeV 
(off resonance). The other supersymmetric masses are as specified in the text. 

in Fig. El the region in the plane that is consistent with EWB (gray shaded band) 

along with the |rfg" loop | = 1.9 x 10~ 27 e • cm (95 % C.L. limit) curves for various values of the 
first generation slepton masses (assumed degenerate). For a given slepton mass, the region 
in the plane consistent with EDM constraints lies below the dashed line. In the same 

figure, we also plot the |<i 2 ~ loop | = 1.9 x 10~ 27 e-cm curve (solid red line) from two-loop SUSY 
contributions [33j . Several key features emerge from Figs. [7| and |H1 

(i) In the range of /1 and M2 we are considering, \d e \ is dominated by the one- loop contri- 
butions for slepton masses below 1-2 TeV, while the two-loop effects become dominant 
for slepton masses larger than 2-3 TeV. 

(ii) In the case of resonant EWB, which requires the smallest amount of CP violation, the 
electron EDM constraint requires slepton masses to be heavier than 1 TeV. 

(iii) Two-loop contributions to d e imply that EWB cannot occur too far off resonance (see 
Fig. IHJ), even in the limit of very heavy sleptons. 

Additional constraints on higgsino-mediated electroweak baryogenesis do arise from the phe- 
nomenology of indirect dark matter detection in the MSSM, and they are investigated in 
Ref. @. 
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FIG. 8: We plot in the |(^>«|-|/i| plane the region consistent with EWB (gray shaded band), the 
|^i-ioop| _ ^ g x e . cm ^grj y c c l ii m it) curves for various values of the first-generation slepton 
masses (dashed horizontal lines), and the |(ie~ loop | = 1.9 x 10~ 27 e • cm curve. For a given choice of 
mass parameters, the allowed parameter region lies below the EDM curves. The baryon-to-entropy 
ratio is required to be in the range 4.8 x KT 11 <Y B < 9.8 x KT 11 0, and the SUSY parameters 
are as in Table [H 

Before concluding, we emphasize that the constraints implied by Fi gs.pl andEl are specific 
to the MSSM, and that the extensions of the MSSM discussed in Ref. and elsewhere can 
lead to different phenomenological conclusions. In particular, extended Higgs sector models 
with additional scalar degrees of freedom can give rise to a strong, first-order electroweak 
phase transition without requiring a light t^. In this case, resonances in the stop sector may 
enhance the importance of CP-violation associated with the tri-scalar terms {e.g., (Pa), and 
the information provided by the neutron and neutral atom EDM searches would become 
more important than for the MSSM scenario considered here. In addition, we also note that 
there could exist additional, 0(1) corrections to Yb associated with computations of the 
sphaleron rates, bubble profile, and Majorana gaugino transport that we have not addressed 
here. 



V. CONCLUSIONS 

The present work is part of a broader program initiated in ^3] whose goal is to system- 
atically reduce uncertainties in EWB calculations induced by transport phenomena. The 
main new results of this work are: 

• We have calculated the contribution to the quantum Boltzmann equations due to 
decays and inverse decays induced by tri-scalar and Yukawa-type interactions. We 
have performed the calculation in the Closed Time Path formalism to leading non- 
trivial order in the ratios = fi/T, e p = T/u, Ed = v w k cS /u. 

• Specializing to the case of MSSM, we have derived the (inverse) decay rate due to top- 
quark Yukawa interactions, their supersymmetric tri-scalar counterparts, and the soft 
SUSY-breaking tri-scalar interactions proportional to y t . These rates are of 0(a®), 
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and - when not vanishing due to threshold effects - they can be of comparable size or 
larger than the 0(a s ) contributions from scattering processes. 

• We have revisited the fast-Fy approximation [il], which consists in taking the rate Ty 
of Yukawa-induced processes as large compared to all other relevant time-scales. We 
have found previously-unnoticed corrections to the baryon density that enter at linear 
order in the 1 /ry-expansion, whose inclusion shows that this expansion in fact breaks 
down. The approximation is sound in the unbroken phase, where Yukawa processes are, 
indeed, fast on the scale of diffusion processes. But in the broken phase, the rates r^-, 
T# associated with relaxation processes can be as large as or larger than Ty, even after 
including the O(ot s ) contributions to T Y - The interplay of these competing processes is 
quite significant, and a quantitative analysis requires performing a numerical solution 
to the transport equations for realistic, finite values of Ty. For the parameter choices 
we considered, keeping Ty finite can increase Yb by a factor between 20% and 100% 
compared to the Yy — > oo limit. 

• We have updated our previous [l5| analysis of the connection between EDM constraints 
and EWB. Even within present uncertainties, the simultaneous requirement of success- 
ful EWB and consistency with EDM upper limits, poses stringent constraints on the 
size of SUSY CP violating phases and mass spectrum. For example, for any value of 
the CP violating phases, successful baryogenesis and one-loop EDM constraints force 
the slepton masses to be heavier than ~ 1 TeV. Bounds of this type will be sharpened 
by future EDM experiments and can be tested at future collider experiments. 
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APPENDIX A: 1 B AND 1 F IN TERMS OF ONE-DIMENSIONAL INTEGRALS 

Performing the l>l integral in Eq. (J 18)) and the uj^ integral in Eq. (|22j) yields 
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and 
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APPENDIX B: DETAILS OF THE SOURCE DERIVATION 

In this Appendix we give some details of the derivation of the source terms reported 
in Section IIIII We first relate the scalar Higgs chemical potentials fin u d to the Higgsino 
chemical potential fig and then show how to further simplify the final expression by use of 
mass relations among weak doublet partners. 

Recall that the Higgsino chemical potential fig corresponds to the vector charges, ng +t0 = 
\I/ g +t0 '-f ^/ jj+.o, for the Dirac fields introduced in Eq. (|29j). which combine u- and d-type 
Higgsino densities. To determine the scalar Higgs density that is kept in equilibrium with 
the Higgsino vector charge density via gaugino interactions, we examine their interactions in 
the MSSM Lagrangian, written in terms of the Dirac fields ^#+,0, and the four- component 
gaugino fields: 



w+ 



(W + 



1P> 



w 3 



t ' 



B 



B 



(Bl) 



The charged wino field is a Dirac spinor, for which a vector charge density can also be 
defined, while the neutral fields ^ w o B are Majorana spinors, whose vector charge density 
is zero. In terms of these fields, the Higgs-Higgsino-gaugino interactions are: 



C 



HHV 
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V 6+ (H?P L + e^H+P R )* B + * 6 o(H° d *P L - e^H° u P R )* B 
<S> 6+ (-H;*P L + e^H+P R )V w0 + <j>Ho(H° d *P L + e^H° u P R )^ w0 



(B2) 



- 92 

+ h.c. 



*H + {HfP L + e^H u P R )y m + V W+ (H-*P L - e^H+Pjn)* 



c 



The combinations of scalar fields appearing in each term of Eq. (jB2|) tell us which densities 
are kept in equilibrium with the Higgsino densities by fast gaugino interactions. To illustrate, 
consider the second term on the RHS that couples the Hq and B fields to the neutral Higgs 
fields. Using fi B = we see from Eqs. (|2(JI21|) that this term in C HH -y will generate source 
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terms for Ho given by 

Sfjo - 



91 



2 L 
9\ 
2 . 



(jiga - fJ-H°)^F(mjj ,mg,m H o) + (/i^o + ii H o)l F (mjj , m B , m H o) 



where 



fJ, H o 



2 



(B3) 
(B4) 

(B5) 

(B6) 
(B7) 



Similar expressions follow from the other terms in Eq. (jB2|) (assuming the W densities 
vanish). The assumption of "fast" supergauge interactions then leads to: 7 



(B8) 
(B9) 



and 



= Hh + Hh HH d = fih- , (BIO) 

where [i B and \Xh refer to the common chemical potentials for the charged and neutral Higgs 
scalars. Adding together the top and stop sources in Eq. (J32l33j) and using the relations 
(l34IB8IBl(Jj) gives for the Yukawa source for the density T: 

St = -NcVt^T - HQ ~ Hh) x B(A t ;m iR ,m iL ,m H o)+l B (A t ]m iR ,m~ bL ,m H +) 

+1 B Gu; m~ tR , m- tL , m H o ) + 1 B (p; m~ tR , m bh , m H - ) 
+ X F (m k+ , m tL , m- tR ) + X F (m k0 , m bL , m~ tR ) 
+1 F (m tR , m tL , m H o ) + 1 F ijn tR ,m bL ,m H +) 



+1 F (m tR ,m B0 , m- tL ) + Z F (m tR , m g 



+ ,m~ bL 



(Bll) 



-1 B (A t ; m 4fl , m~ tL , m^o ) - T B (A t ; m tfl , m~ bL , m#+ 
- Jp [m tR , m ti , m H o )-l F (m tR , m hh ,m H +] 



' This is tantamount to assuming that g\Tp H is sufficiently large compared to the other transport coeffi- 
cients so that ~( 1 ho ~ 0. We leave for future work an explicit test of this assumption. A comprehensive 
analysis that allows for fijjo ^ (J-jjo should also include the effects of non- vanishing gaugino densities, since 
gauginos play an essential role in this departure from chemical equilibrium. Since the neutral gauginos 
arc Majorana fcrmions and possess no vector current density, such an analysis will require in turn a study 
of the axial vector analog of Eq. © 0] , a task that goes beyond the scope of the present work. 
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m tL 


= m b L 


= m Q 


m h 


= mr 


= m Q 




= fllfjO 

±J -u 


= ™h u 




= m H o 


= m H d 


m H+ 


= m^o 


= \fi\ . 



We can simplify further by noting that the masses of weak doublet partners are the same: 

(B12a) 
(B12b) 
(B12c) 
(B12d) 
(B12e) 

With the notation for the masses introduced here we arrive at our final result of Eq. (|4U|) . 
APPENDIX C: ANALYTIC CORRECTIONS OF 0(1/T Y , SS ) 

In this Appendix we solve the transport equations in powers of 1/Ty ss and show that the 
analytic solutions obtained in the Ty — > oo limit can receive 0(1) corrections for realistic 
choices of all the competing rates (T y ,h,m)- 

The zeroth-order solution in l/Ty^s is obtained by considering the combination of 
Eqs. (jH| that is independent of Ty and T ss . Letting Dh and D q be the diffusion con- 
stants for Higgs and quark superfields, respectively [see Eq. (fI5|)]. letting the densities be a 
function of z = |x + v„,t| (the co-moving distance from the bubble surface along its normal), 
and neglecting small corrections proportional to for simplicity, we obtain: 

D q [2T"(z) + Q"(z)} + D h H"{z) - v w [2T'(z) + Q\z) + H'(z)\ 

~T(z) Q(z 



k T k 



Q 



+ r h{ -,m-s^) (C1) 



where S^ T = + S?^ and f = df/dz. The approximate chemical equilibrium enforced 
by Yukawa and strong sphaleron processes implies that the combinations 

A T Q H 

™r k h 



S, 



2Q T 9(Q + T) 



kn k T k 



B 



tend to zero in the limit T Yss — > oo ; so that we can formally expand in l/F YjSS an d treat for 
bookkeeping purposes 5y ~ 1/Ty and 5 SS ~ 1/T SS . The relations between the Q, T, and H 
densities, up to order l/T YtSS are then: 

Q = , ^ {kB - 9kT \ , (H + k H 5 Y ) - , - kBkQ - , 6 SS (C3) 
V k H (9k T + 9k Q + k B ) K H Y) (9k T + 9k Q + k B ) y ' 

k T (2k B + 9k Q ) . , k B k T 

k H (9k T + 9k Q + jfefl) [n + " H ° Y) (9k T + 9k Q + k B ) ° ss ■ 



Substituting these expressions back into Eq. (jClJ) . we obtain the equation for H: 

v w H'(z) - DH"{z) = -TH(z) + S(z) + 5S(z) (C4) 
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where 



D 



D h A 



D q K 



(K + A) 
A 



k H (K + A) 



s 

K 
A 



A 



-S 



(K+Ar ToT 

9k T k Q + k B k Q 
k H (9k T + 9k Q 



(C5) 



Ak T k B 
k B ) , 



and 
55 



k B (2k T + k Q )(v w S^ 



DA 



km 



Y 



DqSy) 



{A/k H )T- M 8y 



(C6) 



represents a correction to the effective source 5 for the Higgs density H. The functions Sy 
and 5 SS appearing in Eq. (|C6J) are determined by substituting the lowest order solution Ho 
into Eqs. (jH]) and read: 



M*) = - 



1 

rv 
i 



k 



D h H' Q \z) - v w H' (z) 
k Q + 2k T 



H (z) 

k 



11 



r ss k B k B 



[D q H>>{z)-v w H' {z)]. 



(C7a) 
(C7b) 



9fc T + 9A; Q 

Although in the unbroken phase Sy iSS ~ T,hs/^y,ss x iJ ^ in the broken phase, they 
can be sizable, with 8y 3> S ss . 

All previous treatments have neglected the SS term in Eq. (jC4|) and thus find only the 
leading-order solution for H. Then the only 1/Ty^ s effects appear to be the Sy >ss terms in 
Eqs. (|C3|) . However, SS induces 0{l/Yy tSS ) corrections to the density H obtained by solving 
Eq. (jC4|) . which must be substituted back into Eqs. (jC3|) to give the full Q,T densities to 
order l/T YtSS . Using the simplified bubble wall profile as in Ref. (with constant sources 
in the region < z < L w ), the explicit solution to Eq. (jC4|) in the region of unbroken 
electroweak symmetry (z < 0), that drives the weak sphaleron processes, reads 



H<{z) 



dye~^ (S(y) + 5S(y)) + 



<JS<(0) 



1 



v w z/D 



with 



v,„ ± 



4TD 



(C8) 



(C9) 



and <S5<(0) is the value that SS takes at z = approaching from the left. The 0(l/Ty tSS ) 
contributions to H live in the terms containing SS. The largest effect arises from the presence 
of SS(y) inside the integral. The overall size of SS is dominated by the term in Eq. (jC6j) 
proportional to SyT^. Moreover, the typical size of Sy is set by (r H /T Y ) ■ (H /k B ), leading 
to: 



5S 



r r? M L 



^D(T M + T H ) 



(CIO) 
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with rr = A/[kn(K + A)] ~ 0.07 . Using earlier estimates of IV and we find 

SS/S ~ 0.1, indeed a small correction. However, when using r#, and Ty as calculated 
in Ref. and the present work within the CTP framework, we find SS/S ~ 1, thus 
invalidating the assumption of fast Ty rates. 

In conclusion, the large Sy.ss corrections in the broken phase induce large corrections to 
the effective source for the Higgs density, which in turn induce large corrections to Q, T 
themselves. What past treatments have derived correctly are the 1 /Ty, S s corrections to the 
relation between Q,T, and H (that is, Eq. (|C3j) ). but not the corrections to H itself. Yet 
this correction, it turns out, is the biggest piece of all. 
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